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ABSTRACT 

We study the quantal energy spectrum of triangular billiards on a spherical 
surface. Group theory yields analytical results for tiling billiards while the 
generic case is treated numerically. We find that the statistical properties 
of the spectra do not follow the standard random matrix results and their 
peculiar behaviour can be related to the corresponding classical phase 
space structure. 



1 INTRODUCTION 



The classical dynamics in polygonal billiards (in particular equilateral triangles) on a 
spherical surface was studied in to investigate the effect of the positive curvature 
on the classical motion. The structure of phase space for these curved triangular 
billiards turned out to be regular but very complex. As a consequence of the focusing 
mechanism on the sphere the phase space is entirely covered by chains of stable 
islands. Inside each of these islands the motion is elliptic and can be labelled by an 
infinite repeating code according to the sequence of reflections. This situation that 
we could call 'piecewise integrable' is very different from the one corresponding to 
plane polygonal billiards which are either integrable (or pseudointegrable) or ergodic, 
according to the rationality of their inner angles (see, for example, 0). The natural 
question we would like to answer in the present work is how this peculiar dynamics 
is reflected in the statistical properties of the quantum mechanical spectrum. 
There is a lot of evidence, mainly numerical, that integrable dynamics leads to Poisson 
statistics while classically chaotic systems satisfy the statistics given by one of the 
Random Matrix Ensembles (see, e.g., f^). Mixed systems show a combination of 
these extremes that depend on the relative size of chaotic and regular regions |Q . On 
the other hand, a different kind of systems neither regular nor chaotic has been shown 
to conform to intermediate statistics: a simple example is the plane pseudointegrable 
billiard which has been studied together with other models with similar behaviour by 
Bogomolny et al 0. 

The peculiar phase space structure of equilateral triangles on the sphere led us to 
believe that their spectral properties would not conform to any universal statistics. In 
this article we calculate the quantal energy spectrum and the eigenfunctions of these 
systems. For generic triangles the calculation will be performed numerically while 
the spectra of tiling triangles will be derived analytically just by using symmetry 
arguments, as was the case in the classical counterpart. It was shown in |l| that 
tiling triangles are very particular systems that not only are integrable but for which 
only periodic orbits are present. We will see how this peculiarity affects the quantal 
spectrum. In section 2 we present the model and the numerical procedure to calculate 
the spectra. Section 3 is devoted to the spectral properties of tiling triangles. In 
section 4 we derive the level spacing distributions of several generic triangles and 
discuss them in connection with the corresponding classical phase space. Finally 
conclusions are presented in section 5. 
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2 THE MODEL 



Given T an equilateral triangle with inner angle u centered on the north pole of a 
sphere of radius R— 1, we are concerned with the eigenvalue problem 

AipE + EijjE = in T and V'b = on T (1) 

where A is the three-dimensional Laplacian in spherical coordinates 9 and given 
by: 

Id d 1 d'^ 

and the eigenvalues E = '^S, S being the energies, are dimensionless. 
The wave function ipE can be expanded as 

^s = Ecmi^'"(^'Cos^)e'"^'^ (3) 

m 

where m is an integer and F"^ is a solution of the Legendre's equation: 

^ [{l-x^)Fix)] + [j^^+E] Fix) = (4) 

with X = cos6. The functions which have to be regular inside the triangle are 
then proportional to the Legendre's functions of the first kind PJJ^{x), with real u 
such that E — viy + 1) and real argument —l<x<l. These ( m 7^ 0) go to 
at x = 1, while they diverge for x = — 1. 
It will be convenient to define : 



F"'{E,x) 



^^^-^^'^ p-ix) (5) 



\\r(u + m + l) 

These functions are real ior m < i/ + 1 and satisfy the following recursion relation: 



rim+2 



{ly + m + 2){iy - m - 1) 



2(m + 1) 



X 



X^ 



m)(z/ 



m 



\ {p + m + 2){iy -m- 1) 



(6) 



which is used in the numerical procedure to compute the F"^. 

In order to obtain a spectrum free from the trivial degeneracies due to symmetry we 
have to desymmetrize the billiard. The spherical equilateral triangle is invariant under 
the point group C^y and can be cut along the 3 reflection planes into 6 triangular 
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subdomains. Solving the eigenvalue problem in each of these irreducible domains 
with given boundary conditions corresponds to finding a fraction of the spectrum of 
a given symmetry class. In the following we will treat the Dirichlet problem in the 
fundamental triangle. This corresponds to solutions which are odd under refiection 
across all the symmetry planes. Subspectra corresponding to other symmetry classes 
can be obtained by imposing Neumann and mixed boundary conditions. In the case 
we are considering the expansion of Eq. (|) reduces to: 

= E czpF^^i.E, cos 6) sin 3p0 (7) 

p 

that vanishes automatically at = and = ^. 

In order to solve the problem we use an improved version of the 'point matching' 
method introduced in 0. This consists in requiring that the Fourier coefficients of 
the wave function evaluated on the boundary 

In = j dse-^'^'^'/'^^ljE (8) 

should vanish. C is the billiard perimeter and ds the line element. The existence of 
a non trivial solution leads to the condition: 

det[Jn,m\ = 0, with n,m = 1, N (9) 

where 

Jn,m{E) = J rfse-2^™/^F™(E,cos^(s))e^'^'^(^) (10) 
Eq. @ determines the energy levels E. 

The numerical calculation was performed for equilateral triangles with inner angle 
^ < cij < that is fairly large triangles for which no difficulties are expected at the 
corners. Before presenting these numerical results we analyze in the next section the 
particular cases of tiling triangles that can be treated analytically. 

3 TILING TRIANGLES 

As already described in [Q] the classical motion in tiling triangles can be studied by 
following a unique geodesic on the topological surface obtained by sewing together a 
finite number of replica of the original billiard. Since the geodesies are closed curves 
on this compact surface, all orbits are periodic. In other words, every trajectory in 
a tiling triangle is restricted to a one dimensional subspace. The system is over- 
integrable, the compactness of the sphere playing the role of an additional integral. 
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We therefore expect a non generic behavior of the level distribution. The eigenvalues 
and eigenfunctions of tiling triangles can be evaluated by making use of symmetry 
arguments exclusively. Let us consider the equilateral triangle with u; = |, that is, 
a triangle whose vertices coincide with those of a face of an octahedron. Following 
the desymmetrisation scheme of the previous section we cut the triangle into 6 trian- 
gles with Dirichlet boundary conditions and calculate the 'desymmetrised spectrum'. 
Since the domain tesselates the sphere the solutions of Legendre's equation @) have to 
be continuous and one- valued at all points — 1 < cos 6 < 1. This restricts u to integer 
values / and \m\ < I. The seeked eigenfunctions will then be linear combinations of 
spherical harmonics that can be determined by requiring that they should vanish 
on the boundaries. These boundary conditions result in a reduction of the spherical 
symmetry of the problem: the full spherical symmetry group including all proper and 
improper rotations in three dimensions will then be reduced to a point group. For 
this particular triangle it will be the symmetry group of the octahedron . That means 
that in order to determine the eigenvalues spectrum, that is, the allowed values of I 
and their degeneracy A; we should determine how each of the representations of the 
full spherical group Df" may be decomposed into irreducible representations of Oh- 
The number of times the a*'^ irreducible representation of the subgroup Oh is con- 
tained in the representation Df of the spherical group is given by 0]: 

Ca = —T.X*a{R)Xl{R) (11) 

where the sum extends over the elements of the group Oh the order of which is 
Qo,^ = 24. The characters Xa{R) can be extracted from the corresponding character 
table while the characters for the proper rotations in the full rotation group are given 
by : 

sm^ 

Here Re denotes a rotation through 6 about some axis. The character of an improper 
rotation (that is, proper rotation times an inversion) is the same as the character for 
the proper rotation for the + representations (that is, for even values of /) and its 
negative for the — (odd values of /) representation. 

By applying he 24 symmetry operations of Oh it can be easily seen that the eigen- 
functions satisfying Dirichlet conditions in the reduced triangle transform as the one- 
dimensional irreducible representation Fj" of Oh ■ Therefore the degeneracy A; of the 
eigenstate of energy Ei = l{l + 1) will be the coefficient of the representation Ff in 
the decomposition of Df according to Eq.( p]). 
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The first levels of the spectrum are shown in Table 1 with their corresponding de- 
generacies. The eigenfunctions have been obtained by means of diagonalizing the 
projector of the representation Ti of Oh '■ they are the eigenfunctions corresponding 
to eigenvalues equal to 1. Since the representation is one- dimensional the projector 
has been expressed in terms of the characters as: 



Again the sum runs over all the elements of the group, and the D-matrices corre- 
sponding to improper rotations are minus the corresponding to proper rotations. 
It is easy to understand that only odd values of I appear in the spectrum of Table 1 
since the representation Ti appears only in the decomposition of D representations 
with negative (natural) parity. Furthermore, it can be seen that for / > 13 all odd 
values of / are present in the sequence. That means that except for the lowest lying 
levels the level spacing S{1) = £'^+2 — -E^ is proportional to /. In order to keep constant 
the classical mean level spacing S which according to Weyl's law is S = A being the 
area of the reduced triangle, the level degeneracy A/ should also increase linearly with 
The spectrum of the tiling spherical triangle is thus dominated by number-theoretic 
degeneracies. These 'accidental' degeneracies have been studied by Itzykson et al 
1^. They appear in some simple integrable quantum systems like, for instance, the 
harmonic oscillator with rational frequency ratios ( where E^ n = m + n + 1 ) and the 
plane integrable polygonal billiards ( for the equilateral triangle -Em,n = + — mn 
). The harmonic potential was extensively investigated by Berry and Tabor ^ . It 
is an example of overintegrable systems, in which the additional integrals are the 
conmensurability relations. The spacing among adjacent levels is constant while the 
degeneracy increases with the energy . This makes the mean level spacing go to in 
the classical limit and no spacing distribution can be defined. The case of the plane 
equilateral triangle presents more analogy with our case in what concerns the mean 
level spacing S which is well defined and given by Weyl's formula. Also there the 
states are increasingly degenerate and separated by increasingly large gaps. 
In order to study the statistics of the nearest neighbor spacing distribution P{S) we 
now unfold the spectrum by defining the sequence ei = N{Ei) where N{E) is the 
averaged integrated level density given by 



L is the triangle perimeter in the spherical metric. The level spacing distribution in 
the scaled spectrum P{s), with si = e/+2 — Q is plotted in Fig.l (b) for a sequence 
of 7500 levels. The distribution is bimodal with a strong tight peak at s = and 
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(13) 




90h r 




(14) 
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a flat component, giving equal probability to all allowed values of si. Since both 
degeneracies and gaps diverge with the energy this distribution is not defined in the 
classical limit. Asymptotically the peak tends to a 6- function and the flat mode 
extends to infinite values of s with a height going to 0. For comparison in Fig.l (a) 
we show the spacing distribution for a sequence of 20000 energy levels corresponding 
to a desymmetrized plane equilateral triangle with Dirichlet boundary conditions . 
Although the distribution looks more poissonian, as in a generic regular system, the 



situation is analogous to the one in the curved triangle: as shown in |10] degeneracies 



and gaps diverge with energy. But, as pointed out by Berry in |TT[ for the similar 
case of the right plane triangle, this divergence is so slow that this nongeneric level 
structure governed by number-theoretic degeneracies only appears for very high-lying 
states. Summarizing, the spectrum of tiling spherical triangles is non generic and al- 
though a mean level spacing can be defined, no spacing distribution P{s) exists in the 
classical case. This peculiarity is to be attributed to the number-theoretic structure 
of the spectrum . The spectra of the spherical tiling triangles with angles u; = ^ 
and u = Y' can be calculated in the same way, by considering the symmetry group 
of the tetrahedron and of the icosahedron respectively. The sequence of eigenvalues 
corresponding to the desymmetrized triangle with = ^ and Dirichlet boundary 
conditions, that transform as the one-dimensional irreducible representation r2 of 
is presented in Table 1. In this spectrum odd and even values of / are allowed . Except 
for the lowest lying levels all values of / are present, the level spacing S{1) = Ei^i — Ei 
is again proportional to I and the spectrum has the general characteristics stated 
above. 



4 GENERIC TRIANGLES 

The numerical procedure introduced in section 2 was carried out for curved equilateral 
triangles with inner angle u ranging from f to We checked the stability of the 
solutions by varying the number of points on the boundary used to evaluate the 
integral in Eq.(|D and the number of partial waves in the expansion Eq.@. The 
completeness of each spectrum was tested by comparing the calculated cumulative 
level density N{E) with the averaged one given by the Weyl formula (0). This test is 
particularly relevant for triangles close to the tiling ones, for which most of the levels 
are near-degenerate and therefore some of them could be missed in the numerical 
calculation. 

In Fig. 2 we show the energy spectra up to E = 2400 for a family of equilateral 
triangles labeled by u. These are plotted as curves Ey{uj). The spectra at the left 
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and right boundaries correspond to the integrable triangles with uj = ^ and u = ^ 
respectively, which have been evaluated analytically as described in the previous 
section. The levels in both sequences are labelled by an integer / and their multiplicity 
increase with energy. 

All the remaining triangles with arbitrary rational or irrational angle u are non- 
integrable and therefore present no degeneracies in their spectra. Accidental degen- 
eracies are not expected neither: by varying u, we are moving in a one-parameter 



family and, according to Berry |12| and references therein, variations in one param- 
eter are insufficient to produce degeneracies. However, in the region denoted as /// 
we observe a great number of quasi-crossings. Their origin will be explained below. 
The level spacing distributions corresponding to the unfolded spectra of six triangles 
with uj ranging from f to ^ are shown in Fig. 3. Between 1300 and 2500 levels were 
considered for each case, except for case (d) {uj = O.SSQvr) where the huge number 
of quasidegeneracies prevented us to go beyond a few hundreds levels. We observe 
that the distributions depend drastically on the parameter uj and do not exhibit any 
universal behaviour. In order to understand the general features of the distributions 
and their angle dependence we will have to keep in mind the results for tiling triangles 
derived in the previous section and also refer to the classical phase space plot which is 
shown for each case in Fig. 4. As seen in for generic curved triangles the classical 
phase space is covered by chains of elliptic islands of regular motion characterized 
by an infinite repeating code. The multiplicity of these chains increases and their 
size decreases with the period of the code and phase space takes a fractal structure. 
As shown in Fig. 4 (a) and (f), for triangles close to the tiling ones phase space is 
almost dominated by one chain of 3 islands (corresponding to the + + + code ) in 
the case u = ^, and by two chains (corresponding to + + + and H — ) in the case of 
u = The region outside these big domains is entirely covered by 'dust', that is by 
chains of islands corresponding to very long codes and therefore with extremely high 
multiplicity and small area. These islands , although regular and of the same type 
as the dominant ones, cannot be resolved quantum mechanically in the region of the 
spectrum we are analyzing. As we go to triangles far apart from the tiling ones (see 
Fig. 4 (b),(c), (d),(e)) other chains of considerable size appear and we will see that 
some of them can be resolved in the range of energies considered. In order to see how 
this is refiected in the level spacing distributions we now go back to Fig. 3. 
Histogram (a) corresponds to an angle u = O.Sllvr ( region / of Fig. 2) and can 
be seen as a perturbation of the tiling u = ^. The Oh symmetry is broken and 
the fact that degeneracies disappear and gaps between levels do not diverge with 
energy but tend to some finite value refiects in a modification of the distribution 
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of Fig. 1 (a). The distribution Fig. 3 (a) is still bimodal but the delta peak at 
s = is shifted (i.e., there is some level repulsion) and the flat component does 
not extend to large values of s. Case (f) corresponding to a; = 0.6497r (also in 
region I of Fig. 2 ) presents an analogous situation, with a breaking of the 
symmetry. But in this histogram the low lying peak is broader (in fact, there is 
no level repulsion) and the high s - component shows more structure suggesting 
a superposition of two uncorrelated spectra. In both cases the distributions, and 
more specifically the presence or not of level repulsion, can be better understood 
by analyzing the corresponding eigenf unctions. At a qualitative level, we expect 
that eigenstates spanning an integrable region will show rigid and regularly spaced 
spectra. Therefore, in case where a single integrable region dominates the phase space 
the spectrum will show repulsion, since all the eigenfunctions are correlated. On the 
contrary, if several classes of eigenfunctions coexist, the distribution can be thought 
as a superposition of uncorrelated spectra and there will be no level repulsion. In case 

(a) the eigenfunctions are of two types: one class lives in the + + + elliptic islands 
which mostly fill the phase space, the other class extends over the region left outside 
the islands, which is a very small fraction of the space, (see Fig. 4 (a) ). Thus, the 
eigenfunctions of the first class, which are the majority, will be correlated leading the 
level repulsion observed in Fig. 3 (a). In case (f) the eigenfuctions are of three types. 
Apart from the extended class which lives in the region left outside the islands, there 
arc two classes of localized functions: one living in the + + + elliptic islands, the 
other in the H — islands (see Fig. 4 (f) ). Examples of the Husimi representation of 
functions belonging to these three classes are shown in Fig. 5. The two classes of 
localized eigenfunctions living in domains which are comparable in size will be weakly 
correlated : this explains the much weaker level repulsion in distribution of Fig. 3 

The distributions corresponding to a; = 0.5257r and cu — O.GSStt are shown in Fig. 3 

(b) and (d) respectively. Both correspond to values of the parameter in region II of 
Fig. 2 . We are now fairly far away from the tiling triangle: the sphtting among levels 
is big, and the two modes have shifted yielding a single peaked distribution. As we see 
in the corresponding classical phase space plot Fig. 4 (b) for u; — 0.5257r the islands 
corresponding to codes different from the dominating + + + are still small and cannot 
be resolved in the range of energies we are considering. Still, there are two classes of 
eigenstates coexisting: the ones localized in the + + + domain and the extended ones 
which occupy now a considerable fraction of phase space. Therefore two classes of 
weakly correlated eigenfunctions contribute to the spectrum. For cu — 0.6337r three 
classes of functions are present, since the H — islands are large, as shown in Fig. 4 
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(d), and can be resolved. 

Summarizing, we can say that the spectra of generic curved triangles are non uni- 
versal and follow a parameter- dependent intermediate statistics. The analysis of 
the corresponding classical phase space plots suggest that the characteristics of the 
distributions depend on the number of domains that can be quantum mechanically 
resolved in a given energy range. The conjecture is the following. If only one big 
domain is resolved the eigenf unctions in this domain are of the same type and, there- 
fore, correlated. If more than one domain can be resolved, the coupling among states 
living in different domains will be weak, and the spectrum will be a superposition of 
uncorrelated spectra, showing no repulsion. 

Another interesting observation is the existence of a high peak at a small value of s in 
the level spacing distributions corresponding to triangles in region III of Fig. 2. Two 
examples are shown in In Fig. 3 (c) and (d) for angles uj = O.STlvr and u = O.SSQvr 
respectively. This peak is due to the presence of numerous quasi-crossings in this 
region, the lowest of them are seen already in Fig.2. To understand the origin of 
these quasi-degeneracies we have to refer once again to the classical phase space plots 
corresponding to these triangles. As pointed out in [Q] generic non tiling triangles 
might have domains in which all trajectories are periodic (coexisting with the families 
of elliptic islands). For example, in all rational triangles with u > ^ the islands 
corresponding to the H — code are constituted of periodic orbits. On the other hand, 
in triangles with inner angle satisfying 



all orbits of type + + + are periodic with period n. 

In region /// we find several triangles with rational inner angle or inner angle sat- 



global symmetry group, as in the case of tiling triangles, the presence of these periodic 
domains which cover a big fraction of phase space (in particular the corresponding 
to the + + + code) lead to quasi degeneracies among the quantum states localized in 
these domains. This is the origin of the huge peak at s = in Fig. 3 (d). Distribution 
(d) corresponds to an angle close to : now the periodic orbits are of type H — with 
period n = 7. Since the area occupied by the H — islands is smaller, the effect is less 
pronounced: the peak is lower and located at s small but different from 0. 




(15) 
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5 CONCLUSIONS 



In this work we calculated the quantal spectra of generic and tiling equilateral tri- 
angles on a spherical surface. The spectra of the tiling triangles could be obtained 
analytically by using symmetry arguments while for the generic ones we performed a 
numerical calculation based on an improved version of the 'point matching' method 
which allowed us to go up to around 2000 levels. 

We found that the spectra of tiling triangles are dominated by arithmetical degen- 
eracies. This is a property shared by other simple quantum systems, in particular by 
plane square and triangular billiards. In the case of curved triangles the distribution 
of the degeneracies is easily determined: since the level gap increases linearly with 
/, that is, is proportional to a/E and the mean level spacing in the classical limit is 
constant, the degeneracies also increase as a/E. But the question asked in p about 
relating these degeneracies with some 'hidden' symmetry remains open also in this 
case. 

As expected, the spectra of generic triangles of variable inner angle uj show no uni- 
versal behavior in this energy range. We first remarked that for triangles close to the 
tiling ones the level spacing distribution is bimodal : this can be seen as a remnant 
of the level clustering in the tiling systems. More in general, for an interpretation 
of the spectral properties we had to refer to the corresponding classical phase space 
plots. These might be governed by one or more than one island of considerable size, 
where considerable means that it can be quantum mechanically resolved in the energy 
range under consideration. The conclusions seem to be the following. If there is only 
one resolved domain all the eigenfunctions living in this domain will have support 
on tori and the corresponding energy levels will be regularly spaced, giving rise to a 
rigid spectrum. This is the case in triangles with inner angle uj close to | for which 
the + + + islands covers a great area of the phase phase ((a) and (b) in Fig.4). If 
several domains of comparable size coexist different eigenfunctions which may corre- 
spond to levels close in energy are localized far apart from each other and are weakly 
correlated. No level repulsion will be present in the spectrum. 

In order to check these conjectures and to get a quantitative understanding of the 
spectra we have to extend the calculation to higher energies. The higher we go, more 
and more domains in phase space will be resolved and an increasing number of weakly 
correlated families of eigenfunctions will be present, localized in these domains. It 
will then be possible to subdivide the spectrum in subspectra corresponding to each 
type of eigenfunctions (according to their localization in phase space) and study the 
characteristics of each subspectrum. This subdivision is not possible in the present 
work since we cannot reach values of k for which the small islands will be explored 
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with significant statistics. To remedy this situation and check whether in the high 
k regime the superposition of many uncorrelated spectra recovers asymptotically the 
Poisson behaviour, we are at present adapting the scaling method [|l^ to the geometry 
of the sphere. 
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Fig.l - (a) shows the level spacing distribution corresponding to the first 20000 levels 
of a desymmetrized plane equilateral triangle. In (b) the level spacing distribution 
for the first 7500 levels of a desymmetrized curved tiling triangle with a; = |. 
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Fig. 2 Energy spectra for a family of equilateral triangles with f < < ^ • Full 
points indicate the spectra of the tiling triangles. Region I corresponds to the 
neighborhood of tiling triangles, while II and III denote intermediate regions (see 
text). 
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Fig. 8 Semiclassical representation of three types of eigenfunctions in the triangle 
with a; = 0.6497r . 
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Table 1: First levels of the spectrum of tiling triangles with ou — ^ and u — The 
energy of each level is E — l{l + 1) and its degeneracy A; . 
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